J/ip production and elliptic flow parameter V2 at 
LHC energy 



R. Peng 1 23 and C. B. Yang 1 ' 2 

institute of Particle Physics, Hua-Zhong Normal University, Wuhan 430079, China 
2 Key Laboratory of Quark and Lepton (Hua-Zhong Normal University), Ministry of 
Education, China 

3 College of Science, Wuhan University of Science and Technology, Wuhan 430065, 
China 

E-mail: pengru_1204@hotmail.com 

Abstract. We apply the recombination model to study J/ip production and its 
elliptic flow in the region of 10 < pr < 20 GeV/c at LHC energy. We show the 
distribution of J/ip as function of the transverse momentum px and the azimuthal angle 
<f>. If the contribution from the recombination of shower partons from two neighboring 
jets can not be ignored due to the high jet density at LHC, the elliptic flow parameter 
«2 of J/ip is predicted to decrease with px- 
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In the framework of the recombination model, the mechanism of hadronization is 
described as recombinations of thermal-thermal partons, thermal-shower partons and 
shower-shower partonspQ. With the shower parton distributions (SPD) obtained from 
fitting the fragmentation functions (FF) [2j [3], the model has given good agreement 
with the experimental data of the hadron production at relativistic heavy ion collisions 
(RHIC) In Pb+Pb collisions at ^/s]^ = 5.5 TeV at Large Hadron Collider 

(LHC), the initial temperature and energy density higher than those at RHIC offer a 
new domain to study the physics of strongly interacting matter and the properties of 
quark gluon plasma (QGP). The recombination model has also been applied to a large 
Pt range of 10 < pr < 20 GeV/c at LHC and predicted a new phenomenon that the 
proton-to-pion ratio can be as high as 20, much higher than that at RHIC [5]. Charm 
quark, abundantly produced at LHC is one of the main observables which can trace the 
initial phase of the collision and provide information on the possible formation of QGP. 
For this reason, we study J/ip production in central Pb+Pb collisions (0-10% centrality) 
at LHC in this paper. 

In the recombination model the meson production is expressed as the sum of TT 
(thermal-thermal), TS (thermal-shower) and 55 (shower-shower) terms, since there are 
two components of partons, thermal (T) partons and shower (5) partons originated from 
hard partons produced in hard partonic interactions. While for 55 term, we should 
consider the shower partons from one jet and shower partons from two neighboring 
different jets and we denote the one-jet contribution as 55(1) and the two-jet 55(2). 
Thus the momentum spectrum of the meson is described as 

dN M = d{NjZ + Njf + iVgW + iyff (2) ) _ 
d 2 p d 2 p 
It has been demonstrated in Ref. [T] that contribution coming from 55(2) is negligible at 
RHIC because of the small density of jet and the small overlap probability of the different 
jets. So 55(2) term is not considered in the calculation at RHIC, e.g. SS = 55(1). 

In Refs.[31 E] we have calculated J/ip transverse momentum spectra for different 
centralities in Au+Au collisions at ^s^n — 200 GeV at RHIC and the results fit the 
experimental data well. The contributions coming from TT and TS become lower than 
that of SS for p T > 5.8 GeV/c [3]. The details of calculating TT and TS terms are 
given in Refs. (3j EJ [7] where the thermal parton distribution is determined by fitting the 
low-pr data of J/ip. There are two parameters fitted by the experimental data at RHIC, 
fugacity of charm quark 7 C = 0.26 and the flow velocity vt = 0.3c. And we will still use 
these two parameters since there are not information of J/ip production at LHC so far. 
Now we focus our attention on SS component which is the dominant contribution to 
the charmed meson production in the region of 10 < px < 20 GeV/c discussed in this 
paper. 

The basal idea of the recombination model is to describe the fragmentation process 
as the recombination of shower partons created in a jet [2]. The fragmentation process 
of a parton i splitting into a meson M is expressed as the fragmentation function, 
then the recombination of two shower partons from one jet 55(1) is related to FF Df 1 
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[H With the abbreviation p for the meson transverse momentum, the corresponding 
one-jet contribution to the inclusive distribution is 



pdpdcp 



where fragmentation function Df 1 describes the possibility of parton i splitting into a 
meson M. FFs for J / ip and other D mesons can be found in Refs. [21 E] and Refs. [21 El HD] , 
respectively, including Dg^, D J J^ , , D®° , D® s and D® s . Fi(q, 0, c) is the probability 
of a hard parton i with momentum q at azimuthal angle in a heavy-ion collision with 
centrality c. For a hard parton created with the distribution fi(k) at the creation point, 
the initial momentum k changes into q after traversing an absorptive distance £. The 
corresponding distribution is expressed by the momentum degradation factor G(k,q,£) 

Fi(q,0 = J dkkfi{k)G{k,q,0- (3) 

In terms of £, G(k,q,£) can be written as a simple exponential form G(k,q,£) = 
q5(q — ke~^) [TTJ. If P(£, 0, c) is the probability of having a dynamical path length 
£ for a parton directed at 0, the hard parton distribution is obtained after carrying out 
the integration over £ 

Fi(q,(f>,c) = J d£P&4>,c)Ffat). (4) 

In Ref.[ll]. a scaling behavior of P(£, 4>, c) is found for the dependencies on and c for 
pion production. 0, c) can be written as a universal function in terms of a scaling 
variable z = £/£(0, c) 

P(£,0,c) = ^)/£(0,c), (5) 

where £ is the mean dynamical length and the scaling function ip(z) has been given in 
Ref.flT] which will be used in our calculations. In this paper we discuss J /if) production 
in central Pb+Pb collisions at LHC with centrality c = 0.05 (0-10% centrality) and the 
initial transverse momentum spectra fi(k) of hard partons at midrapidity at LHC are 
parametrized in Refs.|12j. 

The contribution which involves shower partons from two different jets is negligible 
at RHIC due to the small overlap probability for two neighboring jets At higher 
energy, such as LHC, because of the high jet density the new component of the two-jet 
contribution SS(2) should be taken into account, which is given by 

dN ss ^ 1 ^ fdqdq' 

i 7i = — 7, / -FiiQ^^jFi'iq ,<f),c)T{q,q) 

pdpdcp p [) p J q q 

I ——Fw{q,q';pi,p 2 )RM{Pi,P2,P)- (6) 
J Pi P2 

In the above expression, F u i is the distribution of shower partons related to the two jets 
and for a meson it is written as 

F„{q,<f ] p x ,p 2 )=Si{^)£${^), (7) 



x 
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and SI (z) is the SPD of shower parton j with momentum fraction z in a jet initiated 
by hard parton i. Another fundamental concept of the recombination model is that 
the recombination process involves quarks and antiquarks, while gluons are converted 
to quark-antiquark pairs in the sea before hadronization [13]. Our view is that a hard 
parton creates a shower of partons that recombine subsequently to form hadrons. Thus 
gluons are not included in the shower partons, e.g., i = u(u),d(d), s(s),c(c), g and 
j = u(u),d(d), s(s),c(c). Rm(pi,P2,p) is the recombination function (RF) for the 
process of the two constituent quarks with momentum pi and P2 to form a meson 
with p. The parameterized results of the SPDs and RF for the charmed mesons have 
been discussed before in Refs.[3J E]. The only function we need to know is the overlap 
function T(q, q') between the two neighboring jets, which reflects the probability for the 
overlap of the two shower partons. It depends on the density of jet, or on the energy 
■s/s]^, the momentum vectors q and q' of the hard partons, and the width of their jet 
cones. Since we have not sufficient information of such dependencies for collisions at 
LHC, the overlap function is approximated by an average quantity T which varies over 
a wide range T = 10~ n with n — 1, 2, 3 and 4. 

The results of the inclusive J/ip distribution at a given azimuthal angle = 57r/24 
for four values of T are shown in FigJTJ In order to take a clear comparison of the 
contribution from each term we illustrate them (with T = 10~ 2 for SS(2)) in Figj2]where 
TT term is too small to be shown in the figure. As T decreases, the corresponding 
probability of the overlap between two neighboring jet cones becomes small and the 
distributions get close automatically, as shown in FigJTJ In the px range shown, the 
TS term is about four orders magnitude lower than SS{1) and SS{2) because of the 
rapidly decreasing thermal parton distribution. The contribution coming from SS{2) 
decreases with p? faster than that from SS(1) and SS(1) exceeds over SS(2) when T 
is smaller than 10 -3 . Thus the difference between the distributions of the four values of 
T decreases with increasing pt- 

In Figj3] we show the azimuthal anisotropy of the transverse momentum 
distribution of J/ip at given p? = 15 GeV/c as a function of the azimuthal angle. 
We can get the elliptic flow v 2 for different values of T for such a <fi dependence through 

^77 = A(p T )[l + 2v 2 ( PT ) co8(20)]. (8) 
pxapTdq) 

Shown in FigJH is the predicted t>2 of J/^ as a function of pt for different V. The t>2 
value increases with V. And the values of v 2 for V from 10~ 3 to 10 _1 are larger than that 
at RHIC calculated in Ref.fT], where t>2 increases monotonically up to px — 4 GeV/c 
and then saturates at 0.04 for 0-10% centrality. If V is small enough (r = 10 -4 ) so 
that the impact of SS{2) can be ignored, the value of t>2 maintains as a constant about 
0.04 which is almost the same as the saturated constant calculated at RHIC |7j. When 
r increases to lO^ 1 the SS{2) term becomes so large that the other three terms are 
negligible and t>2 also keeps as a constant of 0.07. The main reason is that when V is 
large the contribution from SS{2) plays an important role because there are two terms 
contributing to the azimuthal angle dependence in SS recombination. 
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The results of v 2 for T = 10 3 and T = 10 2 decrease with increasing p?. In 
fact the trend that v 2 decreases with p? has been probed in heavy ion collisions at 
RHIC at p T > 3 GeV/c [HITS] (review in Ref.p]). STAR Collaboration presented the 
results of the azimuthal anisotropy parameter v 2 for Au+Au collisions at \/Snn = 200 
GeV extend the measurements up to pt = 12 GeV/c p2] from the reaction plane 
method and four-particle cumulant method. It shows that v 2 is maximum at pt ~ 3 
GeV/c, saturates at 3 < pr < 6 GeV/c and then decreases slowly up to the highest 
momentum measured. Thus the reliable discussion of the trends in the region of pt > 6 
GeV/c requires multi-particle correlation measurements [18J. The eventual decrease of 
v 2 at high p T has been demonstrated as a consequence of finite inelastic jet energy loss 
by the perturbative QCD theory [19J. Another recombination model which describes 
hadronization via quark coalescence has also been applied to explain the elliptic flow of 
mesons and baryons [20] . 

In Fig|5j we calculate v 2 coming from the term of TS, SS(1) and SS(2), 
respectively. The values of v 2 from TS and SS(1) terms are almost constants of 0.041 
and 0.038. While v 2 from SS(2) increases with p? slightly and at the same time it is 
much larger than those from TS and SS(1). It's a reasonable result. If the meson is 
formed by the recombination of thermal- shower partons or shower-shower partons from 
the same jet, the process is related to only one jet. While two different jets should be 
considered in the recombination of SS(2). Shower partons from two jets may undergo 
more energy loss than from one jet before they are formed into a meson and the energy 
loss depends on the azimuthal angle <fi of the hard parton in Eq.(jl]). Thus the value of 
v 2 associated to two different jets is larger than the others. Also we have discussed that 
the contribution coming from SS(2) decreases faster than that from SS(1) when pt 
increases. So when the contributions to the flow from the terms of SS(1) and SS(2) are 
comparative (T = 10~ 3 , 10~ 2 ), the dependence of the distribution of J/ip will decrease 
with increasing pt, as shown in FigJH 

Finally we predict the production of the other two charmed mesons D° and D s and 
the dependence on the transverse momentum and the azimuthal angle is illustrated in 
FigHwith T = 10~ 2 . 

In summary, we predict the distribution of J/ip and its elliptic flow parameter v 2 at 
LHC energy in the framework of the recombination model. Because of the high energy 
density the recombination of shower partons from two neighboring jets may play an 
important role in the meson production and its influence to v 2 is discussed. And v 2 of 
J/ijj is predicted to decrease with p T . 
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Figure 1. J/tp distribution as functions of the transverse momentum pr for the four 
different values of T. 
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Figure 2. The comparison of three terms for J/ip with T = 10 




Figure 3. The azimuthal anisotropy of Jip distribution for the four different values 

of r. 
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Figure 4. The predicted Vi of J/ip for the four different values of T. 
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Figure 5. w 2 from TS, SS(1) and SS(2) terms, respectively. 
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